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1 Introduction 

The Density Matrix Renormalization Group (DMRG) is an electronic structure 
method that has recently been applied to ab-initio quantum chemistry. The method 
originated in the condensed matter community with the pioneering work of White 
IlllH. Although the earliest quantum chemistry implementations are only a few 
years old, the DMRG has already been used to solve many problems that would 
have been intractable with any other method, and especially, multireference prob- 
lems with very large active spaces. For example, we have used the DMRG to study 
systems ranging from molecular potential energy curves [3] O, to excited states 
of large conjugated polymers ||5] |6), to metal-insulator type transitions in hydro- 
gen chains (7j. In each case, we have obtained accuracies close to the (estimated) 
exact Complete Active Space Configuration Interaction (CAS CI) or Complete Ac- 
tive Space Self-Consistent-Field (CASSCF) result, for active spaces well outside 
the range of traditional algorithms e.g. 100 active electrons in 100 active orbitals 
0. Unlike a traditional CAS (where the active space wavefunction is obtained in a 
brute-force Full Configuration Interaction expansion) the DMRG utilises a compact 
wavefunction ansatz. However, this ansatz is very flexible, is well-suited to nondy- 
namic correlation, and in the cases of long molecules, provides a near optimal, local 
description of multireference correlations. 
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Historically, the DMRG was not originally formulated from a wavefunction 
ansatz perspective, but rather in the Renormalisation Group (RG) language of Wil- 
son's Numerical RG IH[9j[l]|2), f rom which it is descended. The original quantum 
chemical implementations of the DMRG were also described from an RG point 
of view (e.g. |[T0] [TT] [121 [T3l [14]| ) . Although the mathematical form of the DMRG 
ansatz has been known for some time lfl5l[T6l[T7l[T8l . only in recent years has it been 
realised that the wavefunction view of the DMRG provides a more convenient and 
in many cases more powerful paradigm, and this has led to fundamental advances 
in the DMRG method itself lfT9ll20ll2T1l22ll23ll24ll25l I26ll27ll28l |29ll30ll3Tl 1321171. 

The current article provides an expository introduction to the DMRG in quantum 
chemistry from the wavefunction point of view. This is complementary to earlier 
articles that use the RG based formulation and the first-time reader will benefit from 
reading such articles alongside the current one. It is not our intention to provide a 
comprehensive review of the DMRG method even within the restricted domain of 
quantum chemistry. Thus we do not pretend to survey the literature except to say 
at the start that the field of quantum chemical DMRG has developed through the 
work of White et al.flTOj [33l [331, Mitrushenkov et al.flj]|35]|36l, our contributions 

dHEZlHlEHlTJillil.the W ° rk ° f Legeza ' Hess et al '02] EU gO] S2, the work 
of Reiher et al. l42l l43l l44l [141 . and most recently the work of Zgid and Nooijen 
B31 . Also related, but too numerous to cite in full here, are the developments with 
semi-empirical Hamiltonians; some representative early works are those in l46ll47l 
SUSSES ED |52). In addition > we mention again that the DMRG has its origins in 
the condensed matter community and thus excellent sources of information which 
provide this perspective are the recent reviews of Schollwock [32| and Hallberg 

E3ED. 

The structure of our article is as follows. We begin by introducing the underlying 
DMRG ansatz and examining some of its special properties in sections [2] and [3] In 
sections @] and E we explain the connection between the wavefunction ansatz, and 
the original Renormalisation Group language within which the DMRG is usually 
described. In section [6] we describe how the structure of the DMRG wavefunction 
allows the efficient evaluation of Hamiltonian matrix elements. Finally, we finish 
with some brief thoughts and conclusions in section|7j 



2 Motivation for the DMRG Ansatz 

The primary challenge in quantum chemistry is to find a good approximation to the 
electronic wavefunction of a quantum state. We can express any A^-electron wave- 
function in a complete basis of Slater determinants, through the Full Configuration 
Interaction (FCI) expansion, 
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\y ) = £ f™-"*|„ in2 „3...„,), (1) 

{« i } = {|0),|l a ),|l^>,|2^}}, (2) 
J>i=tf- (3) 



Here \n\ . . is the occupation number representation of the Slater determinant 
where n,- is the occupation of site (i.e. orbital) i. The total number of orbitals is k and 
N is the total number of electrons. 

The dimension of the coefficient tensor f in the above expansion is 4 , which 
is intractable for values of k much larger than 10. Therefore, we would like to find 
an ansatz where f is expressed more compactly. In particular, we would want such 
an ansatz to require only a polynomial amount of information as a function of the 
number of orbitals in the system, k. 

A very simple ansatz would be to approximate the high-dimensional coefficient 
tensor f by a tensor product of vectors iff 1 . . . \jf k , which we shall call site functions, 

f w y/ 1 ®^ 2 ®^ 3 ...®!/. (4) 

Using the notation y/" 1 to denote the nth element of l/r, i.e. y/" ] = y/,', we can also 
write 

Note that each site function y is not an orbital but rather a vector of length 4, and 
\j/ ni , \j/" 2 represent elements of the different vectors y/ 1 ,^/ 2 . This ansatz contains 
only 4k parameters and is certainly tractable. However, it is also not, in general, 
very accurate. So, let us try to improve the ansatz by increasing the flexibility of the 
site functions y/. We can introduce additional auxiliary indices, i.e. 

W np ^Wj- (6) 

The new indices i,i' are auxiliary in the sense that they do not appear in the final 
coefficient tensor f and must be contracted over in some fashion. The simplest 
arrangement is to contract the indices sequentially from one \]f site function to the 
next, i.e. 

f™-"^ E <<i 2 v$ 3 "-v£ 1 - (7) 

hhh—k-i 

For simplicity, we will assume that the dimensions of all auxiliary indices are chosen 
to be the same, and we shall call this dimension M. Then each site function y/ 
is a 3-tensor of dimension 4xMxM, and the total number of parameters in the 
wavefunction ansatz is AM 2 k. 

This is, in essence, the DMRG ansatz for M states. (More precisely, it is 
the ansatz used in the one-site DMRG algorithm, as explained later). Note that 
by increasing the dimension M, we can make the approximation arbitrarily ex- 
act. Because (for given n\ ...n^) the contraction in Eq. (0 is a series of matrix 
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products, this ansatz is referred to in the literature as the Matrix Product State 
Q5l[l![l3[IE|53[l9l|25l|54l|55). Combining the site functions explicitly with 
the Slater determinants we have 

hh'3--ik-l 

Before continuing, let us first establish some notation. The site functions y/ in 
Eq. © are 3-tensors. However, the notation of linear algebra is designed primar- 
ily for vectors (1 -tensors) and matrices (2-tensors). Naturally, any 3-tensor can be 
considered as an array of matrices, so long as we specify which two indices are the 
matrix indices and which is the 3rd (array) index. When viewing the site function 
as an array of matrices, we will write the 3rd (array) index on the top. Thus in this 
notation, we have 

Matrix : [l/A] (dimension M xM) 
Elements : w" p (9) 

and the DMRG wavefunction ([8]) is written as 

l%MRG>= E W n ]W n W']---[W nk ]\nin 2 m...n k ) (10) 
mn 2 n 3 ...n k 

(Note that the first and last site functions [y/" 1 ], [y"*] have dimensions 1 x M and 
Mxl respectively). 

Alternatively, we can view a 3-tensor as a single matrix if we group two indices 
together to make a compound index. This view will be useful when discussing the 
renormalised basis and canonical representations of the DMRG wavefunction in 
sections |4] and [5] Depending on the context, we will either group the n index with 
the left or the right auxiliary indices, giving 

Matrix : [y/ p ] (dimension AM x M) 
Elements : w p . ., 
or Matrix : [\)f p } (dimension M x AM) 
Elements : y? ? (11) 

Note that the superscript p here denotes the pih site function in the DMRG ansatz 
(O, not any particular element of the site function. 



3 Properties of the DMRG ansatz 



Let us now examine some properties of the DMRG ansatz. 
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Variational: Since we have an explicit wavefunction, the expectation value of 
the energy provides a variational upper bound to the true energy and in practice 
DMRG energies are evaluated in this way. As M is increased, the DMRG energy 
converges from above to the exact energy. 

Multireference: There is no division into occupied and virtual orbitals, all or- 
bitals appear on an equal footing in the ansatz (0. In particular, the Hartree-Fock 
reference has no special significance here. For this reason, we expect (and ob- 
serve) the ansatz to be very well-balanced for describing nondynamic correlation 
in multireference problems (see e.g. I37l l4lf71). Conversely, the ansatz is ineffi- 
cient for describing dynamic correlation, since this benefits from knowledge of 
the occupied and virtual spaces. 

Size-consistency: The DMRG ansatz is size-consistent within a localised ba- 
sis. Consider a system AB composed of two spatially separated, non-interacting 
subsystems A and B. Associate localised orbitals 1 ... a with subsystem A and 
a + 1 . . .a + b with subsystem B. Then, the DMRG wavefunction for AB fac- 
torises into a product of DMRG wavefunctions for A and B. First expand the 
DMRG wavefunction 



«rg>= e ^^..<;_ li xtr--Kr,-ti»i---»^+^) 

"1— n a +b 
h—ia+b-1 

= E( E vC-vCJ"! •••«*> 

ia n\...n a 
k—ia-i 

- e ^xX;t>°+i •••««+»>)• < 12 > 

n a+l-- n a+b 

'a+l—'a+b-i 

Then note that we can write a separable wavefunction |f AB ) = |f A )|f B ) for- 
mally as |f AB ) = Zi=i l^ A >l^ B > and thus we can take the dimension of index 
i a which couples systems A and B above to be 1, giving 

Mrg>= E <•.•¥£>•••««> E CCl"«+""^ 

n\...n a n a+l ...n a+b 
h—ia-l ia+l—ia+b 

— I ^dmrg) I ^dmrg) ■ (13) 

Compactness and efficiency of the ansatz: The number of variational parameters 
in the DMRG ansatz is 0(M 2 k). How large do we need M to be to achieve a 
good accuracy? If we choose, for a given index i p , M = 1, then the wavefunction 
factorises into a simple product of contributions from the spaces {n\...tip\ and 
{n p+ \ . . .nit}. Increasing M then introduces additional correlations or entangle- 
ment between the wavefunction components in the two spaces. The M required 
for a given accuracy thus depends on the correlations in the specific state of the 
molecule. However, we have seen in our applications that for appropriate prob- 
lems, even modest M = (9(100 — 1000) can allow us to obtain very good accuracy 
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and to solve problems that are insoluble with other techniques. Of course, having 
a small number of variational parameters does not guarantee that an ansatz can be 
manipulated efficiently. (Witness the difficulty in evaluating the variational en- 
ergy corresponding to a Coupled Cluster wavefunction!) As we shall see in Sec. 
[6] the product structure of the DMRG ansatz enables matrix elements to be eval- 
uated without ever reconstructing the DMRG coefficients in the full Slater deter- 
minant expansion, thus bypassing the exponential complexity. (Although one can 
do so if one wishes, e.g. for the purposes of analysing the DMRG wavefunction, 
as in |[T4l ). Finally, we note that the DMRG incorporates correlations between 
orbital spaces in a sequential manner, i.e. the first set of auxiliary indices i\ en- 
tangles spaces {«i} and {n2--.nj l }, ii entangles spaces \n\ni\ and {n^.-.n^ 
and so on. For this reason, the DMRG ansatz performs best if strongly-correlated 
orbitals are placed next to each other in the ansatz [12, 39l l42ll34l . 
5. A local multireference ansatz for long molecules: The DMRG wavefunction is 
particularly well-suited to long molecules where it can be viewed as a naturally 
local multireference ansatz. In long molecules (i.e. those where one of the di- 
mensions is much larger than the other two) with a finite electronic correlation 
length, we can divide the molecule at any point along the backbone and expect 
the degree of entanglement between the two resulting subsystems to be indepen- 
dent of the point of division and the length of the chain. Thus, for such problems, 
the M required for a given accuracy is independent of the length of the system 
and the number of variational parameters in the DMRG wavefunction is simply 
const x 0(k), as should be in a local ansatz. However, unlike in other local cor- 
relation approaches the DMRG provides a local multireference ansatz. It is this 
local nature even in the presence of strong nondynamic correlations which has 
allowed us to solve very large active space multireference correlation problems 
in long molecules E] . 

In problems which are large in two or three dimensions, the degree of entangle- 
ment between two subsystems grows exponentially with the length of the border, 
and thus the preceding considerations no longer apply. We might then ask, can we 
modify the DMRG ansatz to obtain a naturally local multireference description 
for large systems with arbitrary dimensionality? Recently, this has been shown to 
be possible. Consider, for example, two rows of atoms (each with one localised 
orbital) arranged as in Fig. Q] The first sub-figure illustrates the sequential cou- 
pling between orbital spaces that is contained in the DMRG wavefunction, which 
is inefficient at describing correlations between atoms in different rows. In the 
second sub-figure, however, we have added additional auxiliary indices to cou- 
ple the site functions both along the rows as well as along the columns in a 
non-sequential manner. This is the basis for the so-called Pair-Entangled Product 
State wavefunctions which present one of the most promising new developments 
in this area Il2ni22ll23l . 



An Introduction to the DMRG Ansatz in Quantum Chemistry 



7 




(fh) (jb (jh) 
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(a) DMRG 




(b) PEPS 



Fig. 1 Density Matrix Renormalisation Group and Pair Entangled Product State wavefunctions 
for two rows of atoms. Note in the DMRG ansatz, the site functions are coupled sequentially, 
which prevents the efficient description of correlations between the rows. However, in the PEPS 
ansatz, addition indices are added to the site functions (e.g. — ► y/J^ ^ ) whose coupling directly 
captures the inter-row correlations. 

4 The Renormalized Basis 

As we have discussed above, the auxiliary indices of the site functions introduce 
couplings between the orbital spaces in the DMRG ansatz. In addition, they can also 
be provided with a direct physical interpretation. Just as the index n, is associated 
with the Fock space of orbital ;, so can we also associate a set of renormalised many- 
body spaces with the auxiliary indices of each site function y/. This provides the 
Renormalisation Group (RG) interpretation of the DMRG wavefunction. Consider, 
for example, the first set of auxiliary indices i\. We first perform the summation 
in the DMRG wavefunction expression over n\, which couples l^" 1 with the set of 

states {\n\)} = {|0), |l a ), \2"P)}. This formally defines a space {/i} with basis 
functions |z'i) 

|ii> = lYCI«i> < 14 > 
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or more succinctly 

{h} = f 1 -{n x }. (15) 

Of course, the transformation of the {n{\ orbital Fock space by the y/ 1 site function 
is trivial. (Indeed, if, as is usual, we do not allow \]f to mix states with different 
particle numbers or spin, we would simply have = \n\) for all 4 states). However, 
things are more interesting, when we consider the spaces associated with later sets 
of auxiliary indices. For example, repeating the above exercise for 22 

l'2)=E< 1 <l 2 |«i» 2 ) (16) 

n\H2 
h 

=E<y'i»2), (17) 

"2 

{h} = V 2 ■ {hn 2 } = V 2 ■ V 1 ■ {n\n 2 }. (18) 
In general for the space {i p } and the associated basis \i p ), we write 

{i P } = V ■ {ip-Wp} 

= <jfP-<jfP- 1 ---^-{n 1 n 2 ...n p }, (19) 

"p 

ip-l 

= E ^v^UC^" 2 ---^- (20) 
«i •■■«/> 

Note that the matrix representation of \jf p is simply the matrix form of the site func- 
tion [\j/ p ] described in Eq. (fTTTl . i.e. 

(ip\r\i P -mp) = Yf piip _ im (21) 

and thus we can also write Eq. (|20T > as 

\i P )= E [V n W 1 ]---W , "- l W]\nxn 2 ...n p ). (22) 

n\...n p 

Now the dimension of the i p index and {i p } space is fixed to be at most M in 
the original ansatz ([H). Thus, the action of \[f p ■ ■ ■ iff 1 is a projective transformation 
from the full many-body space down into a renormalised many-body space of M 
basis states, where each basis state \i p ) is expressed as a linear combination of many 
product functions \n\ . . .n p ) with coefficients given by Eq. (1201 . The renormalised 
spaces have a recursive structure: {i p } is obtained from {i p -i} which is obtained 
from {/p-2} and so on. 

The construction of one renormalised space from the previous one may be con- 
sidered to proceed in two stages. To construct the space {i p }, first we form the 
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(a) One site ansatz 
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Fig. 2 Block diagrams for the one and two site DMRG ansaetze. 




(23) 



(24) 



The first step is called "blocking" and the second step "decimation" in the tradi- 
tional language of the Renormalisation Group, and therein lies the basic connection 
between the DMRG ansatz and its RG interpretation. It is common to represent 
these blocking and decimation steps in the pictorial fashion shown in Fig. [2] 



5 The Canonical Representation and Sweep algorithm 

The DMRG wavefunction is invariant to a class of transformations of the site func- 
tions y/, since the associated nested many-body spaces {i p } are themselves invariant 
with respect to transformations within each space. The original DMRG algorithm, 
which was formulated in the language of orthogonal projective transformations fol- 
lowing Wilson's Numerical Renormalisation Group, in fact corresponds to particu- 
lar choices of representation of the site functions within the above invariant class. 
We shall call such representations "canonical representations". All existing DMRG 
implementations in quantum chemistry work with canonical representations of the 
DMRG wavefunction. In addition, the use of canonical representations is closely 
linked with the density matrix interpretation of the DMRG and also with the DMRG 
sweep algorithm, which provides a natural algorithm to optimise the DMRG wave- 
function. 

Associated with each DMRG wavefunction f there are k canonical representa- 
tions, one for each site. At site p, the canonical representation is written as 
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|«P>= £ [L" 1 ] . . . [L"''- 1 ] [C' 1 '' ][/?"''+ '] . . . [R nk ]\n\ ...n p ...n k ) (25) 

n\...n p ...n k 

= E ^•••^ I ^ I ^i-"^ 1 l»i-"«P"-n*>. (26) 

'l — ip-l Sp—fk- 1 

Here, the site functions to the left of p have been given the symbol L, while those to 
the right have been given the symbol R. The L and R site functions, which are in this 
context usually called transformation matrices, are each orthogonal matrices when 
written in the matrix representation of Eq. dTTb . We interpret the L site functions as 
matrices by grouping the n index with the first auxiliary index, 

{q<p:L\ nV :=L n ^ (27) 

and in this form we have 

[m T [L«] = [\\, (28) 

L Lq inAi" = 8 '''"- w 

In 

For the R site functions, we group the n index with the second auxiliary index 

(q>p:R c i rn :=R'X) (30) 

and in this form we have 

[/? 9 ][fl ? ] r = [l], (31) 

L R *,m R km = S ">- (32) 
rn 

The L and R matrices each define a set of orthogonal projective transformations, 
which give rise, respectively, to two sets of renormalised spaces {/} and {r} asso- 
ciated with the site p representation of the DMRG wavefunction. The {/} spaces, 
{'2} ■ ■ ■ are built up by incorporating the orbitals in the order 1,2 . . .p, 

{q<p):{l q }=L q -{l q -in q } 

= Li-li- 1 -{l^ 2 n q - 1 n q } 

= L«-L«- x ---U-{n l ...n q } (33) 

and the |/} functions form an orthogonal renormalised basis (from the orthogonal 
nature of the [L] transformation matrices) for each {/} space 

\l q )= £ [L n '][L^]...[L n »-'][L"»]| ni n 2 ... ni ), (34) 

n\...n p 

(l q \l q ) = S W . (35) 
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The {r} spaces and \r) basis functions are defined similarly, but now the orbitals 
are incorporated "backwards" in the order k,k— 1 ...p+ 1 

(q>p)- {r q }=^-{n q r q+ i} 

= R q -R q+l -{n q n cl+l r q+2 } 

= & -R q+1 ---R k -{n q ...n k }, (36) 
\r q )= E [R" q W q+1 ] ■ ■ ■ [^-'P^h^+i ...n k ), (37) 

n q ...n k 

{r q \r' q ) = 8^. (38) 

Having defined the renormalised spaces, we now see that the C p site function 
gives the wavefunction coefficients in the product space formed from the renor- 
malised left basis the orbital space {n p }, and the renormalised right basis 
{r P } 

\y) = L C lnr\ l P-inpr P ) (39) 

inr 

where we have used the notation Cf„ r : = Cf! p r . 

inr ip-l r p 

We now consider the DMRG wavefunction expressed in the canonical represen- 
tations of sites other than p. Since the same wavefunction is simply being expressed 
in a different representation, this implies a relationship between the wavefunction 
coefficients C and transformation matrices L,R at different sites. Comparing repre- 
sentations at sites p, p + 1 we see 

|f)= £ [L ni ] . . . [L n p- 1 ][C" p ][/?"'' +1 ][/?''p +2 ] . . . |«i ...n p ...nic) (40) 

n\...n p ...n k 

= £ [L m ] . . . [L"p- 1 ][L n ''][C n P +1 ] [R"p +2 ] . . . [R" k ] \ni...n p .. .n k )- (41) 

n\...n p ...n k 

This implies 

[C n p][R"p+i] = [L n P][C n P+ l ] (42) 

or, switching to the alternative matrix interpretation of Eq. (fTTl i for C p ,C p+l and 
likewise for L'\RP +1 

Y^C/'n.rK.r'n = H L 'ln.l' C l' .n> V ' ^) 

r I' 

From C p , we can determine the quantities in the site p + 1 canonical form that do 
not explicitly appear in the site p canonical form, namely C p+1 ,LP, by the singular 
value decomposition (SVD) of C p , 
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Fig. 3 The DMRG sweep algorithm with the one site ansatz. After all L blocks are constructed 
going from R, the sweep direction is reversed. 

Clr = ^ ln vG V V Vr , (44) 

/' 

cCr-L^K^- (45) 

r> 

The connection through the SVD between the representations at different sites 
leads to the density matrix formulation of the DMRG. Recall that the singular vec- 
tors of a matrix M may be related to the eigenvectors of M T M and MM T . Thus from 
C p , we can define a symmetric object [r p ] = [C P ] T [C P ], i.e. 

r in,l'n' = ^L^fnr^Vn'r- 
r 

r p is none other than the density matrix associated with the left subsystem, or 
"block" of orbitals 1 . . .p, and the left transformation matrix LP is obtained as the 
matrix of M eigenvectors 

E r in.l'n' L rn'.I" = L L,l" °f" ' ( 47 ) 
l'n' 

This corresponds to the traditional density matrix interpretation of the DMRG: to 
obtain the canonical representation at a new site requires a basis change into the 
eigenvectors of the subsystem density matrix. 

The sequential set of transformations from representation to representation along 
the sites also yields a natural optimisation procedure for the DMRG wavefunction 
known as the sweep algorithm. At each site p, we solve the Schrodinger equation 
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in the basis {l p _in p r p } to obtain the coefficient matrix C p , thus (dropping the sub- 
scripts on the basis functions for simplicity) 

(l'n'r'\ti-E\>¥)=0, (48) 
£(Z'nV|/? -E\lnr)Cf nr = 0. (49) 

Inr 

From this coefficient matrix, we obtain the new transformation matrix at site p ± 1 
from the SVD in Eq. (l44l (or equivalently, in the density matrix formulation, from 
the eigenvectors of the density matrix in Eq. d47li). If we move through the sites from 
left to right (p —>■ p + 1) in a sweep, we successively determine new LP matrices, 
while moving from right to left (p — ► p — 1) determines new R p matrices. After 
the sites are traversed in one direction, we traverse in the opposite direction thus 
allowing improvement of all the LP and R p matrices. (Of course, to initialise the 
procedure, requires some starting guess for the LP and R p matrices). This is the 
basic method that is employed to optimise the DMRG energy. 

We usually depict the canonical representation at site p in a block-configuration 
diagram as shown in Fig. |2]consisting of a left block of orbitals 1 ...p — 1, the site 
p and a right block of orbitals p + 1 . . . k. Then, moving from one site to another 
corresponds to moving from block-configuration to block-configuration, sweeping 
from left-to-right and then right-to-left as shown in Fig. [3] 

So far we have always been working within what is known as the one-site DMRG 
algorithm, since, as can be seen from the block diagram in Fig. [21 there is only one 
site between the left and right blocks. However, in earlier formulations of the DMRG 
algorithm it was common to use the so-called two-site algorithm, corresponding to 
the second block configuration in Fig. [2] Here the wavefunction at site p is written 
in the renormalised product space as 

1*0 = L C 'Lr\ l P-^ P n p+l r p+l ) (50) 

Inn' r 

where we see that two complete orbital Fock spaces {n p },{n p+ i} appear in the 
wavefunction expansion. Unlike in the one-site configuration, we can only approxi- 
mately relate the canonical representations of the two-site wavefunctions at different 
sites, and thus there is no single consistent DMRG wavefunction across a two-site 
DMRG sweep, but rather a whole family of DMRG wavefunctions, one at each site. 
Originally, the two-site algorithm was introduced to eliminate some numerical prob- 
lems associated with local minima when optimising the DMRG wavefunction in the 
sweep algorithm Il2l [l2l[37ll , but with the introduction of newer methods which avoid 
such minima l56l . the one-site formulation should now be viewed as preferred. 
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6 Evaluation of matrix elements 

For completeness, we now outline briefly how the DMRG wavefunction allows the 
efficient evaluation of the matrix elements necessary to solve the Schrodinger equa- 
tion in the renormalised product basis (|49| |. We first note that any operator in the 
complete Fock space {n\ } ® {ni} ® {nj} ® • • • <8> {«/t} can be expressed as a sum of 
products of "local" operators that each act individually in the Fock space of a single 
site. For example, the unit operator / in the full Fock space may be considered as a 
single tensor product of local unit operators 

t = P ®P®P®...®I k (51) 

where e.g. (n\ \I y \n\) = 8 nin i . To see how the quantum chemistry Hamiltonian 

H = a/ + Y, v ijkia] a ) a k a l (52) 

ij ijkl 

can be written as a sum of products of local operators, it is sufficient to show that 
the creation and annihilation operators can be expressed in this form. Note that a 
single creation or annihilation operator does not simply act in the Fock space of a 
single orbital, because of the anticommutation relations between operators. Instead, 
we write for aj,a; 

4 =Y\(~) n ' ®P ia ]Pi, (53) 
a i = Y[{-) n i®P i a i P i . (54) 

j<i 

Here the operator Y[j<i(~ formally keeps tracks of the anticommutation, since if 
we consider e.g. a, acting on a determinant, it counts the number of sign changes 
involved in moving orbital i to the front of the orbital string. P,- denotes the projection 
of operator onto the {«, } space alone. 

Given that all operators can be written as a sum of products of local operators, 
we now examine how the matrix elements of a single product of local operators are 
obtained. Consider the product 

6 = 6 1 ®6 2 ®...6 k . (55) 

In terms of the product basis {\l Pl n p r p )} of site p, we can write (dropping the sub- 
scripts on the basis functions for simplicity) 

(l'n'r'\d\lnr) = (l'\O l ® . . .®6 p - l \l){n'\d p \n){r'\6 p ®...®6 k \r) (56) 
= {l'\d L \l){n'\d p \n){ r '\d R \r). (57) 

It is sufficient to demonstrate how the matrix elements (l\Oi\l') are calculated as 
those for Or are obtained in a similar manner. From the recursive definitions of the 
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renormalised basis functions |/), \l') in Eqs. (|34] i, ((37) we have 
(l'\0 L \l)= £ [L^][L n %..[L n P-%Ol ,0\ ^...0 p -' , )\L<][L%..[L n '^]. 



n\...n p i 



n\ ...n„ 



(58) 

These multiple transformations may be efficiently organised into groups of two step 
procedures (corresponding to the familiar blocking and decimation steps of the RG). 
Writing 

ol iir (i l \d l \i[)=L>;;ol in ,L';; (5 9) 

the blocking step corresponds to 

0)p ®Ol„, -» (O l 2 ), , , (60) 

while the decimation corresponds to the transformation into the renormalised basis 
({hn 2 } -> {fa}) 

E ^ 2 (^0 2 ) /in2/ ; 4 4-(^0 2 ) /2Z ,. (61) 

Each such transformation has the cost of a matrix multiplication i.e (9(M 3 ), and be- 
cause of the recursive structure of the transformations, the complete matrix element 
(1'\Ol\1) may be efficiently evaluated as a sequence of matrix products with a total 
cost 0{M 3 k). 

For complicated operators such as the quantum chemical Hamiltonian which 
consist of sums over many products of operators, it is clear that there are intermedi- 
ates which can be reused and saved. For example, the matrix elements of a^a^agaio 
and a{fl2fl4fl5 both involve as an intermediate the renormalised representation of 
a[a 2 , which may be stored and reused. In practice, therefore, the optimal imple- 
mentation of the DMRG algorithm in quantum chemistry requires an efficient or- 
ganisation of intermediates and this is primarily where most of the complexity may 
be be found. The interested reader is referred to the literature for further details e.g. 

nainiiniisiiHisi. 



7 Conclusions 



In this article we have attempted to introduce the Density Matrix Renormalisation 
Group (DMRG) primarily from the view that it provides quantum chemistry with 
a new kind of wavefunction ansatz. Consequently, we can analyse and manipulate 
the ansatz in the way to which we are accustomed in quantum chemistry. By ex- 
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amining its structure we arrive at an intuitive understanding of the strengths of the 
DMRG method e.g. in multireference problems, or in long molecules, where it is a 
naturally local multireference approach. A striking feature of the DMRG ansatz as 
compared to other quantum chemical wavefunctions is the recursive structure. This 
is the connection between the DMRG wavefunction and the traditional language of 
the Renormalisation Group, and provides the central mechanism behind the efficient 
evaluation of matrix elements in the method. 

Traditionally quantum chemistry has understood electronic structure in terms of 
the many-electron wavefunction. We hope that by thinking about the DMRG in this 
language, it will not only become more accessible, but new possibilities will arise 
for cross-fertilisation between quantum chemical techniques and the Density Matrix 
Renormalisation Group. 
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